Strong magnetic coupling between an electronic spin qubit and a mechanical resonator 
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We describe a technique that enables a strong, coherent coupling between a single electronic 
spin qubit associated with a nitrogen-vacancy impurity in diamond and the quantized motion of 
a magnetized nano- mechanical resonator tip. This coupling is achieved via careful preparation of 
dressed spin states which are highly sensitive to the motion of the resonator but insensitive to 
perturbations from the nuclear spin bath. In combination with optical pumping techniques, the 
coherent exchange between spin and motional excitations enables ground state cooling and the 
controlled generation of arbitrary quantum superpositions of resonator states. Optical spin readout 
techniques provide a general measurement toolbox for the resonator with quantum limited precision. 

PACS numbers: 07.10.Cm, 42.50.Pq, 71.55.-i 



Techniques for cooling and quantum manipulation of 
motional states of nano-mechanical resonators are now 
actively explored. Work in this field is motivated by ideas 
from quantum information science [J B, testing quan- 
tum mechanics for macroscopic objectsjjjjjij and poten- 
tial applications in nano-scale sensing [5J, |6( . Approaches 
based on mechanical resonators coupled to optical cavi- 
ties 0] , superconducting devices [1, Q or cold atoms [13] 
are presently being investigated in experiments. 

In this paper we describe a technique that enables a 
coherent coupling between the quantized motion of a me- 
chanical resonator and an isolated spin qubit. Specifi- 
cally, we focus on the electronic spin associated with a 
nitrogen-vacancy (NV) impurity in diamond [l~0 ] which 
can be optically polarized and detected, and exhibits 
excellent coherence properties even at room tempera- 
12j . Since its precession frequency depends on ex- 



ture 

ternal magnetic fields via the Zeeman effect, single spins 
can be used as magnetic sensors operating at nanometer 
scales [HEl. 

The essential idea of the present work can be under- 
stood by considering a prototype system shown in Fig.[T] 
Here a single spin is used to sense the motion of the mag- 
netized resonator tip, that is separated from the spin 
by an average distance h and oscillates at frequency uj r . 
These oscillations produce a time- varying magnetic field 
that causes Zeeman shifts of the spin qubit. Specifically, 
the shift corresponding to a single quantum of motion is 
A = g s ^iBG m ao, where g s ~ 2, /lib is the Bohr magneton, 
G m the magnetic field gradient and an = v/fi/2mav the 
amplitude of zero-point fluctuations for a resonator of 
mass to. For realistic conditions, h w 25 nm, w r /27r « 5 
MHz, a w 5 x 1CT 13 m and G rn w 10 7 T/m we find that 
A/2-7T can approach 100 kHz. Such a large shift can be 
easily measured within a fraction of a millisecond by de- 
tecting the electronic spin state More importantly, 
the coupling constant A can considerably exceed both 
the electronic spin coherence time (T2 ~ 1 ms) and the 
intrinsic damping rate, k = u> r /Q, of high-Q mechani- 
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FIG. 1: A magnetic tip attached to the end of a nano- 
mechanical resonator of dimensions (I, w, i) is positioned at a 
distance h ~ 25 nm above a single NV center, thereby creat- 
ing a strong coupling between the electronic spin of the defect 
center and the motion of the resonator. Microwave and laser 
fields are used to manipulate and measure the spin states. 



cal resonators. In this regime, the spin becomes strongly 
coupled to mechanical motion in direct analogy to strong 
coupling of cavity quantum electrodynamics (QED). 

Before proceeding we note that coupling of mechanical 
motion to several types of matter qubits, ranging from 
Cooper pair boxes to trapped atoms, has been considered 
previously [3, 1(1 IH. The distinguishing feature of the 
present approach is that working at nano-scale dimen- 
sions allows us to combine a well isolated spin qubit with 
a large interaction strength, thus enabling the strong 
coupling regime. In what follows we describe how this 
regime can be accessed in the presence of fast dephasing 
(T 2 * ~ 1 /lis) of the electronic spin due to interactions with 
the nuclear spin bath by using an appropriate dressed 
spin basis. We then show how it can be applied to cool- 
ing and quantum manipulation of mechanical motion. 

In the setup shown in Fig. [T] the nano-mechanical res- 
onator is described by the Hamiltonian H r — ftw r a)a 
with uj r the frequency of the fundamental bending mode 



and 



,t 



the corresponding annihilation and creation 



operators. Motion of the magnetic tip produces a field 
|.Btip| — G m z, which is proportional to the position op- 
erator z = an (a + a') and results in a Hamiltonian 



He 



Hnv + tko r a) a + h\(a + cr)S z 



(1) 



Here H^y describes the dynamics of the driven electronic 
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FIG. 2: (a) Level diagram of the driven NV center in the 
electronic ground state, (b) Dressed spin basis states for 
Q ~ |A| and A < 0. (c) Energies of dressed states in the 
presence of external perturbations of the form H nuc = h5„S z 
for fi/|A| = 0.1 (dashed) and fi/|A| = 1 (solid). 

spin and S z is the z-component of the spin operator which 
we here assume to be aligned with the NV symmetry axis. 

The electronic ground state of the NV center is an S = 
1 spin triplet and we label states by \m s ), m s — 0, ±1. 
Spin states with different values of \m s \ are separated by 
a zero field splitting of u>q/2it ~ 2.88 GHz. For moder- 
ate applied magnetic fields [l6j], \hbB\ -C hu> , static and 
low-frequency components of magnetic fields cause Zee- 
man shifts of states | ± 1) while microwave (mw) fields of 
appropriate polarization drive Rabi oscillations between 
|0) and the exited states | ± 1) as shown in Fig. [2]^ a). In 
a frame rotating with the frequencies of the mw fields, 

H NV = + (|0><**l + KX Q 1) . ( 2 ) 

<=±i 1 

where A± and Q± denote the detunings and the Rabi 
frequencies of the two microwave transitions. For sim- 
plicity we restrict the following discussion to symmetric 
conditions, i.e. Aj = A, = il (e.g. using a single 
mw-field and B —* 0). Hamiltonian |(3J) then couples 
the state |0) to a 'bright' superposition of excited states 
| b) = (| — 1) + \ + l))/y/2, while the 'dark' superposition 
\d) = (| — 1) — | + l))/v2 remains decoupled. The result- 
ing eigenbasis of Hnv is therefore given by \d) and the 
two dressed states \g) = cos(9)\0) — s'm(9)\b) and |e) = 
cos(0)\b) + sin(0)|O), with tan(20) = -V20/A. Cor- 
responding eigen- frequencies are u>d — ~ A and oj e / g = 
(-A±VA 2 + 2fi 2 )/2. We will mainly focus on the regime 
A < where \g) is the lowest energy state (Fig. |2Jb)). 

To achieve a resonant coupling, values for O and |A| 
can now be adjusted such that transition frequencies be- 
tween dressed states, e.g. LUd g = u>d — Ug, become compa- 
rable with the oscillator frequency u> r . Rewriting Hamil- 
tonian ([I]) in terms of \g), \d) and |e) we obtain 

H s = huj r a)a + hw eg \e)(e\ + hw dg \d)(d\ 

+ h(X g \g)(d\ +A e |d)(e| + H.c.)(a + a t ), 

where X g = — Asin(0) and A e = Acos((9). Under res- 
onance conditions, u> r « \oJ g d\ Hamiltonian ([3]) reduces 
to the well-known Jaynes Cummings model (JCM), and 



describes coherent oscillations between states \n)\g) and 
|n — where \n) denotes a phonon number state. 

To observe the coherent dynamics associated with this 
model the vacuum Rabi frequency X g must be compared 
to the motional decoherence rate 7 r and random shifts 
of the transition frequency, Au>d g , due to hyperfine in- 
teractions with the nuclear spin bath. The condition 
A 5 S> 7 r , Auid g then defines the strong coupling regime. 
While in principle 7 r = k at zero temperature, we iden- 
tify below 7 r = ksT/HQ as the relevant decoherence rate 
for experimentally accessible temperatures T. Interac- 
tions with the nuclear spin bath are characterized by a 
typical strength S n ~ l/r 2 * ~ 1 MHz which exceeds A. 

To show how the strong coupling regime can be 
achieved, we now study the driven NV center in the pres- 
ence of hyperfine interactions, H nuc — g s ^BB ntZ {t)S z , 
where B n (t) is the effective magnetic field associated 
with the nuclear spin bath. B n (t) is quasi-static on the 
timescales of interest [l6[ but has a random magnitude 
on the order \B n<z \ = hS n / '(<? s /xb). I n Fig. [^c) we plot 
dressed state energies of H' NV — Hj^v + H nVLC as a func- 
tion of 5 n . For f2 — > we recover the linear Zeeman 
shift for the bare spin states | ± 1). In this regime large 
random shifts of u>d g would prevent resonant interactions 
between the spin and the resonator mode. However, for 
Q, ~ |A| where the operator S z has only off-diagonal ma- 
trix elements in the dressed state basis, perturbations 
are suppressed for f2 3> \6 n \. In particular, the quadratic 
shift of the transition frequency LUd g is given by 

AJV = S 2 Jco gd (1 - tan- 2 (tf) + sin 2 (0)) . (4) 

We find that in the present three level configuration not 
only we can eliminate linear shifts of the transition fre- 
quency, but at a particular value of 9 w 6q ~ 0.227T even 
quadratic corrections vanish. At this 'sweet spot' the re- 
maining frequency shift is Aw^ w 0.6 x S^/uj gd . In other 
words, operation at the sweet spot allows us to suppress 
the effect of -ff n uc by several orders of magnitude and 
treat remaining corrections as a small perturbation. 

Example. As an example we consider a Si-cantilever [j| 
of dimensions (I, w, t) = (3, 0.05, 0.05) [im with a funda- 
mental frequency of u} r /2ir « 7 MHz and ao ~ 5 x 10~ 13 
m. A magnetic tip of size ^100 nm and homogeneous 
magnetization M » 2.3 x 10 6 T//io @] produces a mag- 
netic gradient of G m ~ 7.8x 10 6 T/m at a distance h ~ 25 
nm away from the tip and results in a coupling strength 
A/2tt w 115 kHz. For a temperature of T = 100 mK 
and Q values of 10 5 the heating rate is r y r /2iT w 20 kHz. 
Operating close to the sweet spot 9 « 9q and assuming 
S n ~ 1 MHz we obtain (X g ,j r , Aiv gd ) = 2tt x (70,20,2) 
kHz. Hence, this combination enables us to access the 
strong coupling regime of Hamiltonian ([3]). Higher val- 
ues of Q ~ 10 6 and/or a reduction of the system dimen- 
sions to h ~ 10 nm would allow further improvements. 
The strong coupling regime can then be reached even at 
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temperatures up to a few K. 

Applications. Hamiltonian ([3]) allows a coherent trans- 
fer of quantum states between the spin and the resonator 
mode, which in combination with optical pumping and 
readout techniques for spin states [l2| provides the basic 
ingredients for the generation and detection of various 
non-classical states of the mechanical resonator. Here we 
discuss in more detail an optical cooling scheme to pre- 
pare the resonator close to the quantum ground state and 
a general strategy for the generation of arbitrary super- 
positions of resonator states. 

Cooling and state preparation techniques rely on a con- 
trolled dissipation of energy which in the present setting 
can be achieved by optical pumping of spin states as 
shown in Fig. [3] a). A laser excites spin states | ± 1) 
into higher electronic levels from where they decay with 
a rate Ti back to the same spin state, or with a rate Tq to 
state |0). Projected on the electronic ground state we can 
characterize the pumping process by a tunable pumping 
rate T op (t) = ft2(f)r /(ri +T ) 2 between states | ± 1) 
and |0) and the branching ratio a = Ti/To. While off 
resonant excitations at room temperature yield a ~ 1, 
the ideal limit a — ► can be reached using resonant 
excitations of appropriately chosen transitions at lower 
temperatures [17| . Including mechanical dissipation of 
the resonator mode the evolution of the system density 
operator p(t) is described by the master equation 



a) electronically 
excited states 



p(t) =i[p{t),H s ] + k ((N th +l)V[a\ + N th V[a)}) p 

+ r op (t) (2?[|o)(i|] + o2?[|*>(i|]) P . 



(5) 



i=±i 



Here T>[c]p := (2cp& — &cp ~ pc^c) /2 and Nth — 
[exp(huj r /kBT) — is the thermal equilibrium occu- 
pation number for a support temperature T. 

To remove thermal excitations we first study cw- 
cooling of the resonator mode. Assuming T op (t) = T op ^> 
A we eliminate the fast dynamics of the spin degrees of 
freedom and study the effective evolution of the mean 
occupation number (n)(t) = Tr{p(t)a^a}. The resulting 
equation is oftheform (n) = —W((n) — (n)o), with a total 
cooling rate W — W op + k and a final occupation number 
(n)o = (nN t h + A£ p )/W. Here the optical cooling and 
heating rates, W op — S(tu r ) — S(— u) r ) and A+ p = S(—u> r ), 
are determined by the fluctuation spectrum 



S(uj) = 2A 2 Re / rir (S z (t)S z ) e i 



(G) 



which describes the ability of the spin to absorb (ui > 0) 
or emit (cu < 0) a phonon of frequency uj. To achieve 
ground state cooling W op ~ 0(A 2 /T op ) must exceed A+, 
and the thermal heating rate "f r = nN t h, where j r ~ 
ksT/HQ for relevant temperatures fc^T 3> fkv r . 

The spectrum S(uj) is plotted in Fig. [3J for the side- 
band resolved regime A < T op -C f2, |A|,av where indi- 
vidual resonances can be assigned to transitions in the 
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FIG. 3: (a) Optical pumping of spin states, (b) Level dia- 
gram of Hamiltonian Hs, Eq. J3). Wavy lines indicate optical 
pumping processes into and out of state \d) for a = (solid) 
and a > (solid and dashed), (c) Excitation spectrum S(ui), 
Eq. (0 for LUdg = w r , uj e d — 0.53 ui r , T op /uj r — 0.01 and a = 1. 



level diagram show in Fig. [3jb). Under resonance con- 
ditions, uj r = ujdg, cooling is dominated by transitions 
l n )lff) ~~ * \ n ~ corresponding to the peak in the 

spectrum at w « oj r . This cooling process is partially 
compensated by heating transitions which can occur for 
non-zero populations p ee , pdd of excited states |e) and 
\d). While p ee ~ 1/2 under strong driving conditions, 
transitions | n) j e) — ► |n + l)|d) are detuned from reso- 
nance by \oJdg ~ w e d | = |A| and do not significantly con- 
tribute to heating. The remaining resonant heating pro- 
cess, — * |n+l)|<7), is proportional to the occupation 
of the dark state, which is populated only by optical de- 
phasing processes. Independent of we obtain pdd oc a 
such that in the limit of ideal optical pumping a — > 0, the 
dark state \d) remains unoccupied, thus enabling ground 
state cooling with a strongly driven spin. 

We derive analytic expressions for S(u>) using the quan- 
tum regression theorem. To ensure resonance condi- 
tions we choose values for A < and which fulfill 
n 2 = uj r / cos 2 (9) and study cooling as a func- 



tion of the remaining free parameter 9 E [0, 7r/4]. For 

(A 2 /r op )W where 



T op <C io r we then obtain W op 



w = 



8cos 4 (6») sin 2 (9) 



[(l+a)+sin 2 (6»)][l+cos 2 (26») + 3cvsin 2 (26i)/4] 



(7) 



This function is maximized for 9 « 0.2 n and for 
we obtain an optimized damping rate of W op — 



a -> 
0.8 x 

A 2 /r op . Numerical values for W op (9) in the presence 
of hyperfine interactions are shown in Fig. HJa). Here 
we find that close to 9 « 9q not only we optimize W op , 
but cooling is also most insensitive to perturbations. For 
given W op 3> n the final occupation number is 



(n)n ~ -tan 2 (6») + — 

W0 2 W QHW C 



(8) 



op 



By choosing r op ~ A and a ~ we find that conditions 
for ground state cooling coincide with the strong coupling 
regime. Therefore, a single electronic spin can be used 
to optically cool the resonator into the quantum regime 
starting from initial temperatures T ~ 0.1 — 1 K. For 
a^0 the apparent intrinsic limit for the present cooling 
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FIG. 4: (a) Optical cooling rate W op in the presence of a 
perturbation H nuc = HSnS z for a = and T op /uj r = 0.03. 
(b) Preparation of a superposition state \ip) as discussed in 
the text where T = Tr{\ip){ip\p(t)} and p x = Tr{\i)(i\p(t)}. 
The sequence consists of 6 cooling cycles of duration T c ~ 
7r/4A 3 + 2.5/r op followed by a state preparation stage with 
T p ~ 7r/(4A 9 ) + 7r/(v8A s ). Gray bars indicate short 7r/2- 
rotations between states |g) and \d). The parameters used for 
this plot are 7r/A ff = 0.01, \ g /uo r = 0.01 and T op /u>r = 0.05. 



scheme, (no) ~ ot/2, can be overcome by employing a 
pulsed pumping strategy as discussed below. 

Once prepared near its ground state, the resonator 
state can be completely controlled via Hamiltonian ([3]). 
The key mechanism that enables such a control involves 
the unitary evolution under the resonant JC Hamilto- 
nian ([3]) when oj r — uJdg- Specifically, for a time t = 
■k /2X g this evolution maps arbitrary spin states onto su- 
perpositions of motional states with zero and one phonon. 
This procedure can be generalized for the generation of 
arbitrary states of the form hM = J^^Lo c «.l n )ls) as P r0 ~ 
posed by Law and Eberly [l9(. The basic idea is that 
we can construct a unitary transformation U such that 
U\ip) = |0). Specifically, can be coherently mapped to 
a state [?/>') with the maximal phonon number reduced by 
one after a free evolution Ujc{t) = cxp(—iHsT/H) fol- 
lowed by a unitary rotation U x (f3) — exp(— i(f3\d)(g\ + 
f3*\g)(d\)), which can be implemented by an additional 
external microwave pulse. For appropriately chosen pa- 
rameters r and (3 this combination removes first all pop- 
ulation from state \M)\g) and successively from state 
\M — l,d). By iterating this procedure one can step by 
step construct a unitary evolution U which maps \tp) 
onto the ground state |0, g). Then, by starting from 
|0, g) the inverse operation U~ 1 will generate the tar- 
get state = J7 _1 |0, p). For a given M the state 
can be generated within a time < M/2X g and a fi- 
delity T ~ (1 — MTrj r /2Xg) x po, where po is the initial 
occupation of state |0, g). 

As an example of this procedure, we consider the gen- 
eration of the state = (|0) — \2))\g)/y/2 starting from 
a pre-cooled thermal state with (n) = 0.5. To prepare the 
initial state \0,g) with high fidelity, independent of a, we 
use a pulsed pumping scheme. First, with f2(i) = the 
spin is optically pumped into the bare spin state |0). In a 
second step ft(t) is turned on adiabatically such that the 



spin is prepared in state \g) while \d) and |e) remain unoc- 
cupied. Finally, for a time tint = 7r /(4A 9 ) the system un- 
dergoes an oscillation between states \n)\g) and \n—l)\d). 
The repetition of this pulse sequence successively removes 
motional excitations and prepares the resonator in the 
state |0, g) with a probability po ~ 1 — (7T7 r /4A ff ). Next, 
a sequence of two mw pulses and two partial swaps is 
used to prepare the cat-like state (|0) — \2))/y/2. Fig. [4] 
shows the results of a numerical integration of master 
equation ([5]) simulating 6 cooling cycles followed by the 
state preparation sequence. This example demonstrates 
that quantum "engineering" of motional states is possi- 
ble using the present technique. Finally, the mapping 
procedure can be used for spin-mediated readout of the 
mechanical motional states. 

In summary, we have shown that a single electronic 
spin qubit in diamond can be strongly coupled to the 
motion of a nano-mechanical resonator. Such a strong 
coupling enables ground state cooling and quantum-by- 
quantum generation of arbitrary states of the resonator 
mode. Potential applications include the use of non- 
classical motional states for improved AFM-based force 
and magnetic sensing techniques H, as well as for tests 
of fundamental theories 
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